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ABSTRACT 


Spatial  acquisition  strategies  and  technologies  are 
analyzed  for  optical  communication  systems.  Theoretical 
analysis  is  followed  by  examples  and  comparisons  to  aid 
system  design.  Examples  are  given  for  state-of-the-art 
technology  appropriate  for  space  applications. 
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INTRODUCTION 


Optical  communications  offer  attractive  advantages  for  intersatellite 
communications  and  their  implementation  has  been  investigated  [1].  Our 
primary  interest  is  in  two-way  links  where  both  satellites  have  a  transmitter 
laser  and  a  receiving  system.  A  crucial  advantage  of  optical  waves  for 
communications  is  their  short  wavelength,  allowing  very  large  bandwldths  and 
very  high  antenna  gains.  This  in  turn  permits  the  use  of  low  power 
transmitter  lasers  without  sacrificing  the  throughput  rate  of  the  channel. 

The  use  of  large  antenna  gains  induces  very  narrow  transmission  beams. 
The  system  performance  therefore  relies  strongly  on  the  correct  pointing  of 
the  transmitter  and  receiver  optics.  The  initial  uncertainty  on  the 
positioning  is  usually  much  larger  than  the  value  required  for  establishing 
the  link.  The  function  of  a  spatial  acquisition  system  is  to  estimate  the 
position  of  the  partner  satellite  to  within  a  desired  uncertainty  bound.  The 
tracking  system  then  upgrades  this  estimate  and  compensates  for  further 
fluctuations  in  time.  Although  the  spatial  acquisition  may  be  coupled  with 
the  acquisition  of  other  parameters  such  as  time  and  frequency,  we  are 
concerned  here  primarily  about  spatial  acquisition  and  will  consider  the 
other  parameters  only  for  as  far  as  they  interact  with  the  spatial 
acquisition. 

The  spatial  acquisition  problem  considered  as  a  whole  is  very  involved 
and  the  performance  may  depend  on  a  large  number  of  parameters.  The  goal  of 
the  analysis  presented  in  this  report  is  to  clarify  the  tradeoffs  between 
various  designs  and  to  serve  as  a  guide  towards  sensible  engineering 


choices.  These  choices  include  the  selection  between  direct  and  heterodyne 
detection  technologies  and  between  parallel,  serial  or  zooming  strategies. 
The  problem  will  be  Idealized,  resulting  in  mathematical  expressions  which 
may  compare  to  the  real  world  figures  only  through  multiplicative  constants, 
but  which  are  suitable  for  comparisons .  The  idealization  of  the  problem 
Includes  ignoring  the  relative  motion,  Ignoring  synchronization  issues, 
discretizing  the  parameters,  and  modeling  as  a  fixed  sample  size  Bayes  test. 
The  rationale  behind  these  hypotheses  is  developed  in  the  preliminaries. 

The  appropriate  technology  must  be  chosen  keeping  also  in  mind  that  the 
spatial  acquisition  unit  is  an  ancillary  package.  It  should  therefore  not 
dominate  the  transmission  system  performance,  weight  and  power  requirement. 
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PRELIMINARIES 


This  chapter  presents  some  notation  and  vocabulary  used  through  this 
report,  addresses  some  hypotheses  and  idealizations  made  in  the  sequel,  and 
outlines  the  structure  of  the  next  chapters. 

1.1  Notations 

The  two  satellites  are  denoted  A  and  B.  These  symbols  appear  as  a  right 
subscript  to  specify  which  satellite  a  parameter  describes.  Right 
superscripts  are  used  to  specify  some  parameters  more  explicitly.  Integer 
arguments  on  the  parameters  designate  the  zooming  stage,  from  0  to  1-1,  when 
performing  a  zooming  acquisition.  A  list  of  the  main  parameters  follows. 

T 

D  Diameter  of  the  transmitter  output  lens  on  B 

B 

R 

D  Diameter  of  the  receiver  input  lens  on  A 

A 

e  The  electronic  charge 

fQ  Central  frequency  of  the  IF  signal 

F  Noise  factor  of  the  photomultlplying  detector 

Hijfc  Hypothesis,  ijk,  meaning  the  signal  is  present  in 

spatial  slot  1,  time  slot  j,  frequency  slot  k 

Hp  Hypothesis  p,  meaning  the  signal  is  present  in 

spatial  slot  p 

hf(t)  Transfer  function  of  the  IF  filter 

Hf(<u)  Frequency  response  of  the  IF  filter 

h  Planck's  constant 
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K„(t,u) 


Number  of  zooming  stages 

Boltzmann's  constant 

Number  of  sensors  of  the  revelver  on  A 

Covariance  function  of  the  process  x(t) 

Distance  between  both  satellites 

Sufficient  statistic  for  the  spatial  slot  1,  time 
slot  j,  frequency  slot  k 

Array  combining  all  the  sufficient  statistics 

■  .  Desired  accuracy  improvement  in  solid 

angle.  Also  the  dimension  of  the  M-ary  detection 
problem  for  the  receiver  on  A 

Mean  value  of  the  sufficient  statistic  when  only 
noise  is  present 

Mean  value  of  the  sufficient  statistic  when  the 
signal  is  present 

U  B 

-  .  Ratio  of  the  uncertainty  solid  angle  to 

the  illumination  solid  angle  from  B.  This  is  also 
the  number  of  time  slots  for  the  receiver  on  A 

■  W/W8.  Maximal  value  of 

A 

Double  sided  spectral  density  of  the  noise 
Power  of  the  transmitter  laser  on  B 


Signal  power  received  by  A,  averaged  over  the  total 
acquisition  time 

Probability  of  an  erroneous  acquisition 
Probability  of  error  in  a  binary  test 
Probability  of  a  correct  acquisition 
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Power  spectral  density  of  the  process  x(t) 

Starting  time  of  the  acquisition 

Time  Halts  for  the  signal  Integration 

Total  spatial  acquisition  time  for  the  receiver  on 
A 

Maximum  frequency  shift  between  the  received  signal 
and  the  local  oscillator  signal 

Bandwidth  of  the  IF  filters 

Short  time  bandwidth  of  the  IF  signal 

Input  Impedance  of  the  photodiode  amplifier 

Quantum  efficiency  of  the  optical  detector 

Optical  wavelength 

Noise  photon  count  rate 

Signal  photon  count  rate 

Background  noise  photon  count  rate  per  spatial  mode 
Optical  frequency 

Solid  angle  representing  the  uncertainty  B  has  on  the 
position  of  A 

Solid  angle  representing  the  resolution  limit  of  the 
receiver  on  A 

Solid  angle  representing  the  desired  uncertainty 
range  on  the  position  of  B.  This  Is  usually  equal  to 

■5 

Solid  angle  of  the  beam  emitted  by  B.  The  power 
density  of  the  beam  is  assumed  uniform  over  ft? 


2 

o  Variance  of  the  sufficient  statistic  when  only  noise 

n  . 

is  present 

2 

a  Variance  of  the  sufficient  statistic  when  the  signal 

0  is  present 

0  U 

A  Half  aperture  of  the  cone  subtending 

Some  abbreviations  and  definitions  are  listed  below. 


APD 

Avalanche  photo  diode 

AWGN 

Additive  White  Gaussian  Noise 

PMT 

Photo  Multiplier  Tube 

PSD 

Power  Spectral  Density 

A  parallel  illumination  consists  of  illuminating  the  whole  uncertainty 
zone  0^  uniformly  from  B. 

A  sequential  illumination  consists  of  scanning  with  a  beam  of  width 

U  T 

£  .  Satellite  A  then  receives  a  signal  burst  of  duration  TA/N^,  and 

only  noise  during  the  rest  of  TA« 

A  parallel  receiver  operation  consists  of  mapping  the  total  uncertainty 
zone  onto  the  receiver  array  in  the  focal  plane  of  the  imaging  system  on  A, 
and  collecting  the  data  from  the  sensors  simultaneously. 

A  sequential  receiver  operation  consists  of  scanning  the  image  of  in 
in  the  focal  plane,  with  the  available  sensors. 

A  zooming  acquisition  consists  of  reducing  ■ 
several  steps.  At  each  step  i,  fi^(i)  is  reduced  to  l^(i+l). 


JJ 


<$<o) 


down  to 


in 
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1.2  Idealizing  the  Problem 


We  discuss  In  this  section  some  Issues  of  concern  for  the  spatial 
acquisition  and  develop  the  rationale  behind  the  adopted  solutions. 

1.2.1  Modeling  as  a  Bayes  M-ary  Detection 

The  acquisition  system  is  faced  with  a  non-linear  estimation  problem, 
since  the  parameter  Is  imbedded  in  the  received  signal  in  a  non-linear  way. 

A  problem  of  this  kind  is  usually  separated  into  two  steps  [2];  a)  the  range 
of  the  parameter  is  divided  into  slots  of  appropriate  width  and  the  receiver 
decides  in  which  of  the  slots  the  parameter  lies;  b)  a  maximum  likelihood 
(ML)  or  minimum  mean  square  error  (MMSE)  estimator  provides  a  closer  estimate 
of  the  parameter.  In  the  present  analysis,  we  assume  that  task  b)  of  the 
estimation  problem  is  assigned  to  the  tracking  system.  The  acquisition  can 
be  therefore  modeled  as  an  M-ary  detection  problem.  Furthermore,  we  consider 
the  simplifying  approximation  that  the  parameters  imbedded  in  the  received 
signal  are  also  discrete,  each  value  corresponding  to  one  of  the  receiver 
slots;  see  Sec.  1.2.5.  We  assign  the  a-priori  probabilities  and  choose  the 
costs  for  minimum  probability  of  error.  The  test  is  therefore  a  Bayes  M-ary 
detection. 


1.2.2  Sample  Size  Choice 

Once  a  spatial  acquisition  procedure  is  initiated,  it  should  always 


succeed.  The  appropriate  statistical  test  may  therefore  be  the  sequential 
test  [3].  In  a  sequential  test,  the  data  is  collected  and  continuously 


tested  until  a  decision  can  be  made  with  a  sufficient  certainty.  The  test 
duration  is  therefore  data  dependent  and  Is  a  random  variable.  We  consider 
in  this  report  fixed  sample  site  tests  however  (Bayes  tests  in  this  case). 

We  evaluate  for  these  tests,  the  mean  probability  of  error  for  fixed  data 
acquisition  times.  The  duration  required  to  achieve  a  desired  mean 
probability  of  error  with  the  fixed  sample  size  test  is  longer  than  the 
average  duration  required  by  a  sequential  test  for  the  same  error 
probability.  We  assume  here  the  ratio  between  both  durations  does  not  depend 
appreciably  on  the  acquisition  strategy  and  on  the  signal  statistics.  The 
fixed  sample  size  Bayes  test  therefore  gives  useful  results  for  comparing  the 
various  strategies. 

1.2.3  Relative  Motion 

The  relative  position  of  the  satellites  may  be  derived  from  ephemerldes 
data,  but  fluctuations  in  the  orbits  induce  uncertainties  of  the  order  of 
10~3  rad  on  both  azimuth  and  elevation.  Diffraction  limited  communication 
beams  may  have  a  divergence  as  low  as  10”5  rad.  This  is  ouch  smaller  than 
the  uncertainty  regions  and  the  desired  improvement  in  accuracy  is  on  the 
order  of  104 ,  when  expressed  as  solid  angles.  The  positions  of  both 
satellites  are  not  fixed  in  time,  but  ephemerldes  give  an  estimate  of  the 
relative  angular  velocity  of  the  satellites.  The  largest  relative  angular 
velocity  is  200  prad/sec  and  may  be  known  to  within  a  fractional  error  of 
10  .  If  ephemerldes  position  and  velocity  data  are  used  by  an  open  loop 
spatial  tracking  system,  a  100  sec  acquisition  will  incur  a  tracking  error  of 
only  2  x  10-6  rad  which  is  smaller  than  the  diffraction  limited  resolution  of 


the  optics  of  either  satellite.  The  relative  notion  is  therefore  neglected 
in  this  analysis,  assuming  ephemerides  velocity  data  Is  available  and  used  by 
an  open  loop  tracking  system  during  the  acquisition. 


1.2.4  Coupling  of  the  Parameters 

Before  the  optical  communication  link  can  be  established,  not  only  the 
optics  mast  be  correctly  pointed,  but  also  the  frequency,  polarization  and 
timing  must  be  acquired  and  tracked  to  within  tight  uncertainty  bounds. 

The  acquisition  of  any  of  the  four  parameters  may  depend  on  the 
knowledge  of  the  other  parameters.  The  optimal  acquisition  scheme  can 
therefore  Involve  a  simultaneous  acquisition  of  all  the  parameters.  Our 
primary  concern  here  is  spatial  acquisition  but  we  will  examine  briefly  the 
coupling  of  the  spatial  parameter  with  the  other  parameters.  The  starting 
time  of  the  acquisition  sequence  and  the  schedule  of  the  acquisition  is 
assumed  to  be  known  a  priori  or  received  through  another  communication  link, 
with  an  accuracy  sufficient  to  assume  that  a  coordinated  two-way  acquisition 
is  possible.  If  direct  detection  is  used  for  the  spatial  acquisition  system, 
polarization  and  frequency  are  Irrelevant  parameters  as  long  as  the  signal 
frequency  is  within  the  bandwidth  of  the  optical  filter  of  the  receiver. 
Spatial  acquisition  is  therefore  decoupled  from  these  parameters  in  direct 
detection.  When  heterodyne  detection  Is  used,  quadrature  polarization 
detectors  may  be  used  to  decouple  the  polarization  parameter.  However,  the 
frequency  mismatch  between  the  received  beam  and  the  local  oscillator  may  be 
large  and  unknown  before  the  acquisition  proceeds,  for  heterodyne  detection. 


Spatial  and  frequency  acquisitions  are  therefore  coupled  in  this  case.  In 
our  analysis,  we  consider  several  signal  models  and  cover  separate  as  well  as 
simultaneous  frequency  and  spatial  acquisitions. 


1.2.5  Idealizing  the  Geometry  and  the  Signals 

We  assume  here  the  illumination  beam  has  a  uniform  power  density  over 
its  opening  and  its  geometry  permits  a  scanning  of  the  uncertainty  zone 
without  any  overlap.  This  assumption  is  unrealistic  but  the  overlap  and 
diffraction  effects  may  be  accounted  for  by  a  single  correction  factor  on  the 
laser  power. 

We  consider  all  bandlimlted  signals  to  have  a  rectangular  power  spectral 
density  and  all  bandpass  filters  to  have  a  rectangular  frequency  response. 
Other  power  spectral  densities  or  filter  responses  are  accommodated  through 
correctly  defining  the  bandwidths  involved. 

1.2.6  Parallel  Receiver,  Sequential  Illumination 


The  acquisition  time  for  a  parallel  illumination  corresponds  to  the 
T 

particular  case  N.  ■  1 
A 

The  sequential  illumination  and  parallel  receiver  strategy  therefore 
provides  the  figures  for  the  other  non-zooming  strategies. 

1.2.7  Average  Illumination  Power 

lnce  the  probability  of  presence  of  A  has  a  uniform  density  over  8^,  the 
beam  illuminating  this  region  from  B  must  investigate  all  the  area  uniformly. 
Averaging  over  the  search  time  TA,  the  transmitter  laser  power  is  therefore 
uniformly  distributed  over  the  uncertainty  zone  C^.  If  0A  is  the  half  angle 

opening  of  the  cone  subtending  then  8^  -  ir02  for  small  6^  and  the  average 

L  2  2 

signal  power  density  received  by  A  is  V^/nQjL  .  The  effective  aperture  area  of 

R  2 

the  receiver  is  ir/4  (DA>  and  collects  an  average  power  given  by 


The  average  signal  photon  count  rate  \s  is  given  by 


1.3 


Report  Outline 


Ue  concentrate  our  analysis  on  parallel  receiver  operation  and  sequential 
Illumination.  From  our  discussion  In  Sec.  1.2.6,  the  results  obtained  for  this 
strategy  are  easily  applied  to  the  other  non- rooming  strategies.  We  discuss 
this  strategy  In  Chapter  2  for  heterodyne  detection,  and  In  Chapter  3  for 
direct  detection. 

When  a  parallel  receiver  structure  Is  used,  the  search  reduces  the 
uncertainty  zone  to  a  solid  angle  flF(l)  -  ^  fF(0).  The  proper  acquisition 
requires  the  number  of  sensors  to  be  M  ■  fl^(O)/ fiP,  which  Is  typically 
104.  If  K  <  M,  the  sensors  may  be  scanned  along  the  image  of  cP(o)  In  the 
focal  plan*',  and  the  results  of  Chapters  2  and  3  may  be  used  when  multiplied  by 
the  factor  M/K.  For  some  detector  technologies,  the  maximum  number  of  sensors 
may  be  on  the  order  of  10.  In  a  typical  situation,  the  parallel  acquisition 
time  Is  multiplied  by  a  factor  of  about  10  and  the  total  acquisition  time  may 
become  much  larger  than  one  minute,  which  is  excessive.  The  strategies 
considered  above  consist  of  a  single  data  collection  part  followed  by  a  single 
decision. 

Another  approach  Is  to  conduct  a  zooming  acquisition.  At  the  lc^  step 
of  the  zooming,  the  initial  uncertainty  zone  fl^(l)  is  mapped  onto  the 
K-sensor  array  In  the  focal  plane  and  is  reduced  to  fl^d+l)  -  1/K  fl^l)  by 
the  decision  made  after  the  partial  search  time  T(i).  If  the  acquisition 
proceeds  In  a  cooperative  two-way  fashion,  each  satellite  illuminates  only  the 
residual  uncertainty  zone  8°d)  at  step  1;  the  signal  to  noise  ratio 
therefore  increases  significantly  at  each  stage.  Zooming  acquisition  Is 
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SPATIAL  ACQUISITION  -  HETERODYNE  DETECTION  -  ONE-STEP  SEARCH 


2.1  Organization 

We  discuss  In  this  chapter  the  implementation  and  the  performance  of  the 
one-way  spatial  acquisition  performed  with  heterodyne  detection,  a  sequential 
illumination  and  a  parallel  receiver  operation.  A  statement  of  the  problem  is 
given  in  Sec.  2.2.  A  mathematical  model  of  the  problem  is  given  in  Sec.  2.3. 
Solutions  are  given  in  Sec.  2.4  for  various  assumptions.  Numerical  results  are 
given  in  Sec.  2.5.  Extensions  to  other  search  strategies  are  addressed  in 
Sec.  2.6. 


2.2  Statement  of  the  Problem 

Let  two  satellites  A  and  B  perform  a  mutual  spatial  acquisition  procedure 

starting  at  the  same  time  tQ.  Both  satellites  have  transmitter  lasers  of 
L  L 

powers  P^,  Pg  respectively,  operating  at  the  same  wavelength  X,  and  have 

T  T 

beamforming  optics  of  diameters  D^,  Dg.  The  on-board  receivers  consist  of 

R  R 

optical  imaging  systems  of  aperture  sizes  D^,  Dg  and  heterodyne  receiver  arrays 
of  K&  and  Kg  elements  respectively. 

The  satellites  are  separated  by  a  distance  L  and  the  a  priori  knowledge 
that  B  has  on  the  relative  position  of  A  is  summarized  by  a  uniform  probability 
of  presence  density  in  the  solid  angle  Conversely,  Bg  is  the  solid  angle 

of  the  zone  which  represents  the  a  priori  uncertainty  that  A  has  on  the 
position  of  B;  see  Fig.  2.1.  From  our  previous  discussion,  we  consider  a 
reference  frame  where  there  is  negligible  relative  angular  motion. 


The  spatial  acquisition  procedure  la  considered  complete  when  the  spatial 
uncertainty  zones  have  been  narrowed  from  their  initial  values  and  £^,  down 
to  the  desired  values,  ^  and  We  consider  in  the  sequel  the  case  where 

the  resoltulon  limit  of  the  receiver  on  B,  We  assume  a  fixed  time 
search,  and  consider  the  cost  criterium  to  be  minimum  probability  of  error. 

The  problem  then  becomes  a  Bayes  hypothesis  testing  problem  because  all 
a  priori  probabilities  and  costs  are  assigned.  We  may  then  derive  optimum 
receivers  and  their  probability  of  error  PE.  Subsequently  we  may  Invert  this 
result  to  obtain  an  expression  for  the  acquisition  time  required  to  achieve  a 
defined  probability  of  correct  acquisition  PAC  ■  1  -  PE. 

We  consider  in  this  chapter,  only  two  simultaneous  but  Independent  one-way 


acquisitions.  Since  the  two  acquisitions  are  equivalent,  we  only  evaluate  the 
acquisition  time  TA,  defined  as  the  data  collection  time  required  for  the 
receiver  on  A  to  locate  the  position  of  B  with  a  probability  PAC.  The  duration 
Ta  depends  essentially  on  the  signal  models,  the  signal  to  noise  ratios,  and 
on  the  dimension  of  the  detection  problem.  The  averaged  signal  level  is  given 


g 

by  P  in  Eq.  (1.1)  and  the  dimension  of  the  problem  is  characterized  by  the 

A 

T 

resolution  gain  M^,  the  number  of  timeslots  and  the  number  of  detectors 

S  T 

K  .  We  find  in  this  chapter,  expressions  for  T,  as  a  function  of  P. ,  M  ,  N  , 

A  A  AAA 

and  the  signal  and  noise  models.  Since  the  same  expressions  hold  for 
Tg,  the  subscripts  will  be  sometimes  omitted. 
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2.3  Mathematical  Model  of  the  Problem 

We  assume  first  that  the  one-way  acquisition  is  performed  with  a  parallel 
processing  in  one  step.  This  means  that  the  number  of  sensors  must  be  no 
less  than  the  desired  resolution  gain  -  ft^/ft^.  If  this  is  not  the  case, 
either  a  sequential,  a  zooming  or  a  combination  of  both  receiver  operations 
must  be  used;  see  Sec.  2.6.  We  assume  in  the  following  discussion  that 
KA  -  M*. 

The  illumination  strategy  consists  of  sequentially  illuminating  from  B, 

the  whole  uncertainty  zone  ft^  with  a  beam  of  solid  angle  ft®  >  The 

A  O  ™  D 

U  TUB 

uncertainty  zone  ft  is  therefore  divided  into  N  ■  ft ./ft.  non-overlapping  cells 
A  A  A  B 

T 

and  the  acquisition  duration  T^  is  divided  into  time  slots.  During  each 
timeslot,  the  beam  illuminates  one  cell  and  then  is  switched  to  the  next  cell 
for  the  duration  of  the  next  timeslot;  see  Fig.  2.2.  The  signal  received 
at  A  is  therefore  a  burst  with  a  duration  ^/N^. 

The  receiver  optics  in  A  map  ftg  onto  the  array  of  M  elements  in  the  image 
plane.  Each  sensor  cell  corresponds  to  the  resolution  limit  of  the  optical 
system  and  has  therefore  the  size  of  the  Airy  disk;  see  Fig.  2.2.  The  incoming 
light  is  mixed  with  a  local  oscillator  to  produce  a  heterodyne  signal.  We 
assume  that  there  is  a  perfect  spatial  match  between  incoming  and  reference 
waves . 

Three  models  for  representing  the  IF  waveform  due  to  the  received  signal 
after  heterodyne  detection  will  be  considered.  All  bandwldths  given  are 


one-sided. 
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8  C°  the  f°cal  Plane  of  the 


le 


.•*»  “•  *>  - , . 


‘mi 


13U4»n 


a)  The  signal  Is  modeled  as  a  zero-mean  narrowband  Gaussian  process 
of  bandwidth  W8  with  a  random  center  frequency  fQ  (constant 
over  time)  which  Is  uniformly  distributed  over  a  frequency 
Interval  W. 

b)  The  signal  Is  a  sine  wave  with  a  frequency  jitter  such  that  the 
instantaneous  frequency  is  uniformly  distributed  over  a  bandwidth 
W8  the  center  frequency  of  which  (fQ  constant  over  time)  Is  in 
turn  uniformly  distributed  over  W  >  W8. 

c)  The  signal  is  a  pure  sinusoid  with  random  phase  and  frequency  fQ 
uniformly  distributed  over  the  bandwidth  W;  see  Fig.  2.3. 

The  noise  at  the  receiver  is  dominated  by  quantum  shot  noise.  It  may  be 
modeled  as  an  additive  white  Gaussian  noise  (AWGN)  of  double-sided  spectral 
height  NQ/2  -  hv/2i>  The  rationale  behind  these  signal  models  is  the 
following.  The  optimum  receiver  for  model  c)  gives  an  upper  bound  on  the 
performance  because  it  is  equivalent  to  assuming  an  instantaneous  frequency 
acquisition  and  an  error-free  frequency  tracking  during  the  spatial 
acquisition.  Model  b)  is  a  closer  representation  of  semiconductor  lasers 
noise,  although  the  actual  power  spectral  density  (PSD)  of  the  frequency  noise 
may  not  be  rectangular.  Model  a)  is  an  approximation  to  model  b)  which  leads 
easily  to  optimum  receiver  designs  and  which  is  completely  characterized  by 
second  order  moments.  This  model  introduces  signal  amplitude  fluctuations 
however,  which  may  not  acccurately  model  the  real  signal  when  the  lasers  are 
very  coherent.  For  all  models,  the  center  frequency  fQ  is  random  due  to  the 
frequency  uncertainties  between  the  transmitter  laser  in  satellite  B  and  the 


local  oscillator  In  satellite  A.  The  mismatch,  due  mainly  to  temperature 
changes  and  doppler  shifts,  may  be  large  so  that  W  >  Ws.  However,  the 
acquisition  time  Is  assumed  to  be  short  enough  so  that  fQ  does  not  change 
during  the  entire  acquisition  process. 

2.4  Solution  of  the  Problem 

2.A.1  Reduction  of  the  Problem 

According  to  the  discussion  in  Chapter  1,  we  consider  only  discrete  values 
of  the  random  parameters.  It  is  therefore  assumed  that  the  incoming  signal 
energy  is  concentrated  into  only  one  of  the  M  sensors,  the  signal  is  a  pulse  of 
duration  T/NT  starting  at  a  discrete  valued  random  time  tQ  +  1  T/N^  and 
the  frequency  spectrum  lies  completely  In  one  out  of  NF  non-overlapping 

frequency  bands.  Under  the  above  mentioned  assumptions,  the  Incoming  wave  may 
T  F 

be  one  of  M  x  N1  x  N  orthogonal  signals,  which  are  equiprobable.  Since 
the  spatial  acquisition  is  not  concerned  with  the  estimation  of  the  frequency 
and  time  of  arrival  of  the  signal,  our  problem  is  a  composite  M-ary  detection 
with  a  set  of  unwanted  parameters  which  may  take  NT  x  values  on  each 
hypothesis. 

The  received  waveform  in  each  spatial,  time  and  frequency  slot  may  be 
characterized  by  a  sufficient  statistic  where  iell.M]  labels  the 

spatial  array  cell,  je(l,N^]  labels  the  timeslot,  ke[l,NF]  labels  the 
frequency  slot.  The  random  variables  are  statistically  independent  and 

have  probability  densities,  conditional  to  the  signal  being  in  spatial,  time 
and  frequency  slots  (p,q,r)  given  by 


(1) 


Piijk|Hpqr<LijkKqr)  "  IPs<Lijk)  "  Pn(Lijk)]  6ip 6jq  +  Pn(Lijk) 

where  p8(0  is  the  density  when  the  signal  is  present  in  slot  (ijk)  and 
pn(«)  is  the  noise  density.  The  received  waveform  is  completely 
characterized  by  the  array  _i  “  {ftjkl*  Pro*>ability  density  of  the 

sufficient  statistics  J.,  averaged  over  the  unwanted  parameter  set,  conditional 
on  the  spatial  hypothesis  Hp  is  given  by: 


NT  /  p  (L  ) 

To  ^  ^  P  (LPq  )  n  Pn^ijk^ 

N  N  q*l  r-1  n  pqr  all  (ijk) 


(2) 


The  optimum  decision  rule  for  minimum  PG,  with  equal  a  priori  probabilities,  is 
the  maximum  likelihood  decision  rule. 


For  high  signal  to  noise  ratios,  the  function  Ps^Lpqr^pn^Lpqr^  has  an 
exponential  dependence  on  Lpqr.  This  implies  that  the  suboptimal  rule: 

psa  r) 

max  [max  max  S,  ]  (4) 

pqr  Pn  pqr 


is  close  to  optimal.  This  rule  is  the  optimal  decision  when  the  costs  are 
assigned  for  the  simultaneous  acquisition  of  position,  frequency  and  time  of 
arrival.  If  the  function  ps(*)/pn(0  is  monotonically  increasing,  the  rule 
in  Eq.  (4)  is  equivalent  to  finding  the  largest  Lpqr: 
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( 


max  [  max  max 
p  q  r 


and  declaring  hypothesis  Hg  If  p  ■  s  Is  the  spatial  cell  corresponding  to  the 
largest  Lpqr.  The  decision  rule  in  (5),  as  opposed  to  (3),  does  not  require 
the  knowledge  of  the  signal  parameters.  The  performance  of  the  suboptimal 
receiver  in  (4)  or  (3)  is  well  approximated  by  the  union  bound  in  the  case  of 
low  PE  and  AWGN.  An  error  is  likely  to  occur  when  one  of  the 
(M-l)  x  NT  x  NF  variables  ijjk  corresponding  to  spatial  cells  which 
receive  only  noise  exceeds  the  variable  which  contains  the  signal  component. 
The  probability  of  error  is  bounded  as: 


PE  <  (M-l)  x  NT  x  /  x  PE(2)  0 

where  PE(2)  i8  the  probability  of  error  in  a  binary  hypothesis  testing  with 
the  same  probability  assignments  and  no  unwanted  parameters.  Our  remaining 
work  in  this  section  consists  of  finding  expressions  for  Jlpqr  as  a  function 
of  the  received  signal  jr(t)  and  computing  the  binary  detection  performance 
PE^2). 


2.4.2. 1  Signal  Model  (a):  Gaussian  Noise  with  a  Flat  PSD  Over  W8 

The  sufficient  statistic  in  this  case  corresponds  to  the  detection  of  a 
zero-mean  Gaussian  process  in  additive  white  Gaussian  noise  and  is  given  in 
[ 4 ] ,  Eq .  (2.31)  as 

Tf 

*iik  "  fT  /  r  (t)  h  (t,u)  r  (u)  dtdu  ;  T.  -  T  -  T/NT  ( 


where  r*(t)  is  the  received  signal  in  the  spatial  cell  i.  The  function 
hi(t,u)  satisfies  the  Integral  equation  ([4],  Eq.  (2.23) 


N 

h^(t,u)  + 


I 

T, 


h,(t,z)  K  (z ,u)dz 
1  s 


Ks(t,u) 


(8) 


for  T  <  t,  u  <  Tf 


where  Kg(t,u)  is  the  signal  covariance  function.  Under  the  SPLOT  (Stationar> 
Process,  Long  Observation  Time)  condition,  Eq.  (8)  is  approximated  by  [4], 


I  f  2 

-i-  /  (r(t)  *  hf(t))Z  dt 

°  Ti 


(9) 


with 

,  l2  S  (<o) 

|Hf(w)|  -  s  (<*))  +  N  /2  (1 

8  O 

Since  our  signal  spectrum  Sa(u)  is  modeled  as  rectangular,  (11)  is  equivalent 
to  within  a  multiplicative  constant  to 

|Hf(u))|  -  1  for  juj  e  lfQ  “  f  .  Vf  1 

0  otherwise. 


The  optimal  receiver  front  end  is  shown  in  Fig.  2.4.  The  probability  density 
of  the  sufficient  statistic  is  shown  in  Appendix  1  to  be  conditionally  Gaussian 
under  the  assumption  W®T/N^  >  1.  When  a  signal  is  present. 
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where  Ps  is  defined  in  Eq.  1.1.  When  only  noise  is  present, 
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The  optimum  receiver  requires  the  implementation  of  NF*  -  W/W8 
frequency  filters  for  each  of  the  M  array  cells.  If  a  suboptimal  receiver  is 
designed  with  NF  <  NF*  filters,  then  expressions  (11)  are  shown  in 
Appendix  1  to  become 
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T  F 
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(1 


2  s 

The  expressions  for  m  and  o  are  obtained  by  setting  P  *  0  in  Eq.  (13).  Wi 

n  n 

this  model,  the  ratio  Ps(L)/Pn(L)  is  not  a  monotonically  increasing 


function  so  that  the  decision  rule  of  choosing  the  largest  &ijk  is  even  more 
suboptimal  than  Eq.  (4)  (see  Appendix  2).  The  "pick  max"  decision  rule  has  a 
probability  of  error  given  in  Appendix  2  for  the  binary  case  by 


Inserting  Eqs.  (13)  and  (14)  in  (6)  gives 


this  approximation  is  accurate  when  UT/N^N^  >  100  and  the  bound  is  tight 
when  the  right  hand  side  is  <  1.  This  last  condition  prevails  when  the*system 
has  good  performance.  The  first  condition  may  not  be  satisfied  when  full 
frequency  processing  is  performed,  »  N^*,  and  the  illumination  beam  has 
the  diffraction  limited  opening  NT  -  M. 


2. 4. 2. 2  Signal  Model  (b) :  FM  Noise  with  Flat  PSD  over  Ws 

Because  the  random  process  is  not  Gaussian,  exact  computations  are  far 
less  trivial  in  this  case.  It  may  however  be  intuitively  argued  that  the 
receiver  front-end  in  Fig.  2.4  must  be  close  to  optimal.  The  probability 
density  of  &ijk  In  the  presence  of  a  signal  in  slot  ijk  is  evaluated  in 


(16) 


i. 

a 
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2N  P  T 
o 


Comparing  Eqs.  (16)  and  (13),  we  notice  that  the  expressions  of  mg  are  the 
2 

same  and  that  In  Eq.  (13)  contains  three  terms.  The  first  term  corresponds 
to  the  noise  fluctuations,  the  second  term  to  the  cross  correlation  between 
signal  and  noise,  and  the  third  term  to  signal  fluctuations.  The  third  term  is 
absent  from  Gq.  (16)  as  may  be  expected  and  both  other  terms  have  exactly  the 
same  expression  in  Eqs.  (13)  and  (16).  Expression  (16)  is  a  general 
formulation  of  the  mean  and  variance  after  energy  detection  in  a  bandwidth 
W/NF  for  a  signal  with  no  amplitude  noise.  The  bound  for  the  probability  of 
erroneous  spatial  acquisition  with  a  signal  which  has  no  amplitude  fluctuations 
is 


PE  <  (M-l)  NV  erfc 
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(17) 


2. 4. 2. 3  Signal  Model  (c):  Random  Phase  Sinusoid 
The  waveform  due  to  the  received  signal  is 


s(t) 
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(18) 


The  sufficient  statistic  is  given  in  [5]  as. 


-  A2  +  R2 
±ik  Cijk  8ijk 


'ijk 


’ijk 


^Tf  -  Ti 


cos  (2xf  t)  r(t)  dt 
o 


sin  (2xf  t)  r(t)  dt 
o 


(19) 


and  has  the  following  conditional  probability  densities 


2  -  f  <l2«k  +  E> 

Ps(Lijk)  “  N  Lijk  e  Xo 
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_  _1_  (l  )^ 

2  No  (  ljk> 

poaijk>  •  TTLijk  * 


W  v) 


(20) 


The  exact  probability  of  error  for  the  receiver  structure  in  Eq.  (5)  is 
computed  in  Appendix  3  and  is  well  approximated  by 


PST 

2N 


PE  <  2WT  (M  -  1)  e 


(21) 


which  is  independent  of  the  illumination  strategy. 
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2.5  Discussion  of  Che  Result 


The  expressions  obtained  in  this  chapter  and  in  the  following  chapters 
apply  for  any  problem  geometry  and  hardware.  The  present  work  was  conducted  as 
a  part  of  the  LASERCOM  project,  and  we  have  therefore  used  the  related 
state-of-the-art  figures  to  generate  the  examples  and  graphs.  The 
communication  lasers  are  solid  state  GaAlAs,  operate  at  wavelengths  of  0.8  pm, 
and  are  limited  in  power  to  approximately  10  mW.  Frequency  stabilized  lasers 
may  attain  short  time  bandwidths  of  about  10  MHz.  A  heterodyne  receiver 
bandwidth  of  1  GHz  is  assumed. 

The  doppler  frequency  shifts  may  be  even  larger  and  must  be  compensated. 
The  maximal  distance  occurs  for  a  communication  between  two  geostationary 
satellites  and  may  be  as  large  as  twice  the  radius  of  the  orbits.  Since  the 
optics  mist  have  diffraction  limited  resolution  and  must  meet  some  weight 
limits,  their  apertures  are  limited  to  about  0.1  m.  We  therefore  generate  all 
the  examples  of  this  report  with  the  following  figures  which  represent  the 
state  of  the  art  for  coherent  detection  hardware. 
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As  a  consequence, 

Nc  -  2.48  10"19J 

P8  5.1  10"  lSV 

M  -  1.05  104 

We  use  the  same  figures  for  comparison  with  direct  detection  In  the  next 
chapters,  although  these  numbers  may  not  reflect  the  state-of-the-art  for 
direct  detection  technology. 

The  time  required  to  perform  the  acquisition  with  PE  *  10-4,  one  way,  is 
computed  according  to  formulas  (15),  (17),  (21)  for  the  three  signal  models. 
Curves  giving  T  as  a  function  of  -  1/NT  are  displayed  in  Fig.  2.5 

for  different  values  of  the  number  of  filters  N**  used.  We  draw  the  following 
comments  on  these  curves.  When  the  signal  model  includes  amplitude 
fluctuations,  the  curves  display  first  a  decrease  of  the  acquisition  time  T 
when  the  beamwidth  S2®  is  decreased.  The  acquisition  time  reaches  a  minimum 
for 


and  then  increases  for  larger  values  of  N  .  This  effect  may  be  attributed  to 
the  increase  in  signal  variance  when  it  is  integrated  over  a  duration 
comparable  to  its  temporal  coherence.  However,  the  Gaussian  approximation  is 
no  longer  valid  in  this  circumstance  and  exact  Chi-square  distribution  figures 


should  be  used  instead 
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Fig.  2.5.  Acquisition  times  for  parallel  receiver  operation 


When  no  amplitude  fluctuations  are  modeled  in  the  signal,  the  function 
T(N^)  is  monotonically  decreasing  and  the  minimum  occurs  at  *  M,  when 
the  illumination  beam  has  the  diffraction  limited  opening.  In  this  particular 
case,  the  search  time  is  5  or  30  ms,  depending  on  whether  or  not  simultaneous 
frequency  acquisition  is  performed.  The  performance  bound  given  by  signal 
model  (c)  is  3.2  ms.  The  acquisition  times  for  PE  -  10“ 4,  for  signal  model  (a) 
or  (b)  are  bounded  by: 


where  a  «  0  for  model  b)  and  a  *  1  for  model  a).  Signal  models  a)  and  b)  give 
markedly  different  results  when  the  third  term  is  dominant  In  (23).  The  actual 
signal  after  the  mixing  process  has  only  slight  amplitude  fluctuations  for  very 
coherent  fields,  so  that  the  result  for  signal  b)  may  be  close  to  the  actual 
performance.  Unless  otherwise  stated,  the  signal  model  b)  will  be  considered 
in  the  following  discussions. 

The  short  search  times  predicted  above  correspond  to  an  order  of  magnitude 
of  100  received  photons  (T  -  5  ms).  However,  these  performances  are  obtained 

|i 

only  if  10  heterodyne  detectors  are  available.  This  hardware  constraint  may 
be  difficult  to  meet  in  practice.  We  will  therefore  derive  in  the  following 
section,  expressions  for  the  performance  when  the  number  of  detectors  K  <  M. 
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2.6  Extension  to  One-Way  Serial  Acquisition 

When  the  number  of  sensors  K&  In  the  receiver  is  smaller  than  the 

desired  resolution  Improvement  M^,  the  cells  in  the  uncertainty  zone  ft” 

may  be  searched  sequentially  by  the  sensors  in  a  scan  of  steps.  The 

number  of  scan  steps  must  be  M^/K^.  The  decision,  made  after  the  complete 
scan  is  performed,  is  identical  to  the  parallel  receiver  decision.  The 
acquisition  time  of  the  parallel  receiver  is  therefore  simply  multiplied  by 
M/K.  Expression  (17)  becomes 


PE  < 


(2 


Figures  2.6  and  2.7  show  curves  of  sequential  receiver  acquisition  times. 
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FM  NOISE  SIGNAL.  SEQUENTIALRECEIVER  OPERATION  HETERODYNE  DETECTION.  NF  =  100,  PE  -  1CT* 


Fig.  2.7.  Acquisition  times  for  serial  receiver  operation:  with  simultaneous 
frequency  acquisition. 


3 


SPATIAL  ACQUISITION-DIRECT  DETECTION  ONE-WAY  SEARCH 


3.1  Organization  of  the  Chapter 

Some  analysis  of  spatial  acquisition  with  direct  detection  has  been  done. 

A  good  survey  of  the  results  may  be  found  in  Gagllardl  [6].  We  therefore  limit 
ourselves  in  this  report  to  list  the  expressions  and  derive  numeric  results 
which  may  be  compared  to  those  given  in  Chapter  2  for  heterodyne  detection. 

This  chapter  addresses  parallel  and  serial  receiver  operations;  zooming 
operation  is  left  to  Chapter  4.  Section  2  handles  the  case  where  only 
background  and  quantum  noise  are  present.  Section  3  considers  the  effect  of 
amplifier  thermal  noise.  Section  4  addresses  the  issue  of  the  excess  noise 
produced  by  photo-detectors  with  multiplication.  Sections  2  to  4  consider  only 
a  parallel  receiver  operation.  Section  5  covers  the  extension  to  serial 
processing  at  the  receiver. 


3.2  Quantum  and  Background  Noise  Limited  Performance 
Let  us  first  consider  parallel  illumination.  The  received  continuous 
signal  power  P  is  given  by  Eq .  1.1.  The  mapping  of  the  uncertainty  zone  on 
the  focal  plane  of  the  receiver  on  A  is  covered  by  M  detectors,  each  covering 
an  area  of  the  size  of  the  Airy  disk  of  the  optical  system.  The  number  of 
detected  photons  obeys  Poisson  statistics,  with  count  rates  Xj,  +  Xs  for  the 
sensor  which  receives  the  signal  and  X^  for  the  other  sensors.  Assuming  that 
we  measure  the  exact  number  of  detected  photons,  the  optimal  decision  rule  is 
to  assign  the  signal  to  the  detector  with  the  highest  count  after  the 


acquisition  duration  T.  The  probability  of  error  PE  is  given  by 
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assuming  that  the  acquisition  is  erroneous  when  a  noise  count  equals  the  signal 
count.  Approximations  to  Poisson  error  probabilities  are  addressed  in 
Appendix  4.  The  Union-Chernov  bound  is  tight  for  large  XbT  and  easy  to 
invert. 


PE  < 


~(  v'X.  + 

(M  -  1)  e 


X 

s 


(2) 


where  X8  *  np8/hv.  Because  each  detector  covers  one  spatial  mode,  we  have 
Xb  “  Xbo  where  Xbo,  the  background  noise  count  rate  per  spatial  mode,  is 
assumed  constant  over  space.  In  a  situation  where  a  geostationary  satellite 
tries  to  acquire  the  position  of  a  low  orbit  satellite,  the  background  field  is 
the  earth.  A  typical  receiver  optical  bandwidth  of  3  A  collects  then  a 
background  count  rate  Xbo  on  the  order  of  10 5  s“ When  looking  toward  the 
high  orbit  satellite,  the  count  rate  Xbo  is  much  larger  if  the  sun  is  in  the 
field  of  view,  and  significantly  lower  otherwise.  However,  the  star  images 
introduce  significant  and  unknown  spatial  fluctuations  in  the  background 
noise.  The  numbers  given  in  Sec.  2.5  and  used  for  the  examples  above  give 
X8  -  20,500  a-1. 

When  sequential  illumination  is  considered,  the  receiver  splits  the 
observation  time  T  into  equal  slots.  The  expected  signal  count  in  the 
signal  slot  is  unchanged  but  the  background  count  is  reduced  by  the  factor 


3R 


PE  <  N  (M  -1)  e 
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where  X],  and  X8  are  given  by 


The  bound  in  (3)  may  be  Inverted  as 
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Curves  of  the  time  required  to  perform  the  acquisition  with  an  error 
probability  smaller  than  10-4  according  to  Eq.  (5)  are  shown  In  Fig.  3.1.  Th» 
slight  upward  slope  In  the  curves  corresponding  to  low  noise  counts  in  an 
artifact  of  the  union  bound.  The  union  bound  overestimates  the  PE  for  a  dark 
background  by  a  factor  of  NTM  as  is  seen  from  the  exact  result  for  Xb  ■  0. 
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For  large  background  counts,  the  Gaussian  approximation 
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is  very  good  as  is  seen  from  the  comparison  of  Figs.  3.1  and  3.2 
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BACKGROUNDS!  UANTUM  NOISE.  UNION-CHERNOVBOUND.  PE  =  10-4 


Klg.  3.1.  Acquisition  times  for  parallel  receiver,  background  and  quantum 
noise  limited  operation,  Chernov  bound. 
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Fig.  3.2.  Acquisition  times  for  parallel  receiver,  background  and  quantum 
noise  limited  operation,  Gaussian  approximation. 


3.3  Background  and  Amplifier  Noise  Limited  Performance 


We  make  here  the  a-posteriorl  verified  assumption  of  large  counts  so  that 
Poisson  statistics  are  well  approximated  by  Gaussian  statistics.  The  received 

2 

charge  (Integrated  current)  is  modeled  as  N(m  ,  on>  when  only  noise  is  present 
2 

and  N(m  ,  o  )  when  the  signal  is  present. 

8  8 

If  we  use  an  amplifier  of  input  impedance  Z,  the  thermal  noise  is  an  AWGH 
of  double-sided  spectral  height  N0/2  ■  kt/Z  where  r  is  the  amplifier  noise 
temperature  in  K.  We  therefore  have: 


m  -  eA.T 
n  b 


a2  ,  J«  X-  +  e2A  T 
n  Z  NT  b 


e(A.  +  A  )T 
b  s 


2  ki  I  ,  2,.  .  . 

oa  -  —  -=•  +  e  (X  +  X.  )T 

s  2  nT  s  b 


Acquisition  times  derived  from  (9)  for  various  values  of  Z  are  shown  in 
Fig.  3.3  for  Ab0  =*  10®.  The  acquisition  time  is  almost  independent  of  Afc0 
in  the  range  Z  -  1  to  10® . 
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BACKGROUND+THERMAL  NOISE,  UNION-GAUSSIAN  APPROXIMATION.  XB0  =  10®.  T  =  300K,  PE  =  lO”4 


Fig.  3.3.  Acquisition  times  for  parallel  receiver,  thermal  noise  limited 
operation. 


3.4  Background  and  Quantum-Limited  with  Excess  Detector  Noise 

From  the  graph  in  Fig.  3.3,  it  is  seen  that  photodiode  detection  is  mostly 

corrupted  by  amplifier  thermal  noise.  Improvement  of  performance  versus 

thermal  noise  is  obtained  by  photomultiplication.  In  a  photomultiplier  tube 

(PMT)  or  avalanche  photo  diode  (APD),  every  detected  photon  produces  G 

electrons.  The  multiplying  factor  G  is  random  however,  and  this  introduces 

what  is  called  excess  noise.  The  statistics  of  G  are  described  by  the  mean 

2  -  2  -  2 

gain  G  and  the  noise  factor  F  *  G  /G  «  1  +  var  (G)/(G)  .  If  the  amplifier 
noise  is  dominated  by  background  and  excess  noise,  then  the  means  and  variances 
in  ( 8)  become 
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As  a  result. 


PE  <  (M  -  1)  N  erfc 


/2  A.  +  A 
L.  b  s 


(10) 


(ID 


The  acquisition  times  derived  from  (11)  are  therefore  just  the  times  derived 
from  (6)  multiplied  by  F.  The  noise  factor  F  is  on  the  order  of  2  for  a  PMT 


and  10  for  an  APD. 


3.5  Sequential  Receiver  Operation 


When  the  number  of  available  sensors  K  is  smaller  than  the  number  of 
spatial  modes  in  the  uncertainty  zone  M,  the  M  spatial  slots  may  be 
investigated  sequentially  by  the  K  sensors  in  a  number  of  steps  given  by  M/K. 
The  decision  rule  is  the  same  as  for  a  parallel  receiver  operation  and  the 
required  acquisition  time  for  a  defined  probability  of  error  is  the  parallel 
receiver  acquisition  time  multiplied  by  M/K.  Figs.  3.4  and  3.5  show  typical 
curves  for  sequential  acquisition. 
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Fig.  3.4.  Acquisition  times  for  serial  receiver,  background  and  quantum  noise 
Limited  operation,  Chernov  bound. 
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ZOOMING  SPATIAL  ACQUISITION 


4. I  Preliminaries 

In  the  preceding  chapters,  the  spatial  acquisition  consisted  of  a  data 
collecting  part  followed  by  a  unique  decision.  We  discuss  in  this  chapter  a 
different  acquisition  strategy  where  the  decision  is  taken  in  several  stages. 
At  each  stage,  a  residual  uncertainty  zone  Is  investigated  and  a  partial 
decision  reduces  the  estimated  range  of  the  parameter.  The  residual  range  is 
then  more  finely  investigated  during  the  next  stage.  The  process  is  repeated 
until  the  desired  accuracy  is  obtained.  We  present  in  Sec.  4.2  the  notations 
used  for  this  chapter.  We  address  in  Sec.  4.3  the  optimal  partition  of  both 
the  total  search  time  and  the  allowed  probability  of  error  at  each  stage.  We 
apply  these  results  with  heterodyne  detection  in  Sec.  4.4  and  direct  detection 
in  Sec.  4.5. 


4.2  Notation  and  Framework 

Each  stage  i,  i  -  0,  1,  ...,  I-lof  the  zooming  consists  of  a  spatial 
acquisition  by  Itself.  The  number  of  required  stages  is  given  by 

I  -  logRM  (1) 

where  M  is  the  desired  resolution  gain  or,  otherwise  stated,  the  number  of 
spatial  modes  in  the  initial  uncertainty  zone  when  diffraction  limited 


resolution  is  desired 


At  each  stage,  the  residual  spatial  uncertainty  zone  Sir(i)  is  mapped 
onto  the  K-sensor  receiver  array  and,  after  collecting  the  received  signal 
during  the  partial  acquisition  time  T(i),  the  receiver  decides  which  sensor 
receives  the  signal.  This  decision  is  taken  with  an  error  probability  PE(i). 

The  noise  statistics  vary  between  the  stages  because  the  number  of  spatial 
modes  covered  by  each  detector,  M/K*  +  *  is  not  constant.  For  heterodyne 
detection,  the  multimode  mixing  efficiency  is  lower  than  for  monomode  operation 
and  the  signal  to  noise  ratio  is  therefore  decreased.  This  phenomenon  is 
modeled  by  multiplying  the  noise  level  by  the  number  of  spatial  modes  covered 
by  the  detector 


N  (i  ) 
o 


-(i+1)  hv 


For  direct  detection,  the  background  noise  count  rate  is  the  product  of  the 
background  noise  count  rate  per  spatial  mode,  which  is  a  constant,  by  the 
number  of  spatial  modes  covered  by  each  detector 

A  (i)  -  M-K~(i+1)A. 

b  bo 


The  remaining  noise  contributions,  i.e. ,  amplifier,  quantum  and  excess  noise 
are  unchanged  for  direct  detection. 

If  one-way  acquisition  is  considered,  then  each  satellite  illuminates  the 
full  initial  uncertainty  zone  ii^(O)  during  the  whole  acquisition.  The 
incident  signal  power  for  parallel  illumination  P8  and  the  number  of  "time 
slots"  Nt  which  is  the  ratio  of  the  illuminated  zone  solid  angle  flU(i)  to 
the  Illuminating  beam  solid  angle  Sl®( i ) ,  are  constant  throughout  the  search 


P8(i)  -  P8(0) 

NT(i)  -  Nt(0) 


(4) 


In  a  cooperative  two-way  acquisition,  only  the  residual  uncertainty  zone 
nu(i)  is  illuminated  at  each  step.  The  parallel  illumination  incident  signal 
power  at  each  stage  is  therefore  increased  as 

Ps(i)  _  ps(0)  Ki  (5) 

Several  strategies  may  be  followed  for  serial  illumination,  to  narrow  the 
illumination  beamwidth  from  the  initial  opening  fiB(0)  ■  Q^(0)N^(0)  down 
to  the  diffraction  limit  flB.  Two  strategies  are  investigated  here.  The 
"wide  beam"  strategy  (A)  corresponds  to  maintaining  the  initial  illumination 
beamwidth  ftB(0)  until  it  fills  the  entire  residual  uncertainty  zone  Jl^(i), 
then  using  parallel  illumination. 


fiB(i)  -  min  lfiB(0),  flU(i)] 
NT  -  max  [NT(0)  if1,  I] 


Strategy  A 


(6) 


The  "narrow  beam"  strategy  (B)  corresponds  to  maintaining  the  same  ratio 
between  the  illumination  beam  ftB(i)  and  the  residual  uncertainty  zone  ft^(i) 
until  the  diffraction  limited  beamwidth  flB  is  attained;  see  Fig.  4.1. 


fiB(i)  -  max  ln®(0)  K_i,  flRJ 
NT(i)  -  min  [NT(0),  MK_i] 


Strategy  B 


(7) 


We  consider  that  the  receiver  operation  is  parallel  at  each  stage.  The 
performance  expressions  derived  in  Chapters  2  and  3  therefore  apply  with  the 
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nb  =  nu  =  nr 

B 


nb  =  ft*1  =  nr 

B 


Fig.  4.1.  Two-way  zooming  illumination  strategies. 


above  modifications.  For  a  two-way  acquisition,  the  light  travel  time  from  A 
to  B  must  be  added  between  the  stages.  This  fixed  additional  time  will  not  be 
considered  in  the  analysis  presented  in  this  chapter. 


A. 3  Optimization  of  the  Zooming  Schedule 

The  problem  to  solve  is  the  minimization  of  the  total  acquisition  time 
1-1 

T  «  l  T(i)  (t 

i«=0 

with  the  constraint  that  the  total  probability  of  error  be  no  more  than  a 
prescribed  value  PE.  This  inequality  constraint  is  clearly  met  with  equality 
at  the  optimum 


PE  -  1  -  H  (1  -  PE(i))  (S 

i-0 

We  approximate  the  above  expression  by 
1-1 

PE  -  l  PE(i)  (If 

i-0 

which  is  a  good  approximation  for  small  PE.  Our  analyses  in  Chapters  2  and  3 
give  the  performance  expressions  for  each  stage  as 

PE(i)  -  fjL(T(i))  (11 

The  optimum  is  attained  when 


af^TCi)) 
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where  X  is  a  Lagrange  multiplier.  We  consider  particular  expressions  of  the 
function  ft(»)  in  Secs.  4  and  S.  We  give  here  two  general  expressions  for 
f ^( • )  and  derive  the  corresponding  optimality  relationships. 


a)  f1(T(i))  -  A(i)  erfc  [B(i)  /?(i)  ]  (] 

where  A(i)  and  B(i)  are  parameters  which  may  differ  for  the  various  stages  1 
but  are  Independent  of  T(l).  For  small  PE,  the  optimal  partition  is 


PE(i) 


l  B(  i  )  ] 


b)  f1(T(i))  -  A(i)  exp  [-B(i)  T(i)] 


(1 


The  optimal  partition  of  PE  is 

red)  -  C 

The  problem  is  solved  by  first  computing  PE(i)  from  Eqs.  (10)  and  (14)  or 
(16).  The  partial  acquisition  times  are  obtained  by  inverting  (13)  or  (15). 


4.4  Heterodyne  Detection 

Expression  (2.17)  gives  the  error  probability  for  heterodyne  detection, 
and  a  signal  with  no  amplitude  fluctuations.  If  we  consider  that  frequency 
acquisition,  if  any,  is  performed  only  at  the  first  stage,  we  have 
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where  1^(1)  -  1  for  i  #  0  and  NT(i),  P8(i),  N0(l)  are  defined  in 
Eqs.  (2),  (5),  (6),  and  (7).  The  expression  of  f^C*)  is  in  the  fora  of 
(13).  The  optimal  partition  of  PE  is  therefore  given  by  (14).  Figure  4.2(a-f) 
shows  the  optimal  partial  error  probabilities  and  acquisition  times,  as  well  as 
the  total  acquisition  time  in  the  case  K  ■  4,  *  1.  Figure  4.3(a-f)  shows 

the  same  graphs  for  «  100.  One-way  illumination,  two-way  "wide  beam" 
illumination  and  two-way  "narrow  beam"  illumination  are  addressed  in  each 
figure.  It  is  seen  from  these  figures  that  the  first  stage  of  the  acquisition 
is  predominant  in  all  cases.  Since  this  stage  is  not  affected  by  the  zooming 
strategy,  the  total  acquisition  times  are  comparable  for  the  three  strategies 
as  is  seen  in  Fig.  4.4. 

4.5  Direct  Detection 

In  the  case  of  direct  detection,  the  expression  for  f^( • )  depends  on  the 
noise  statistics,  which  may  in  turn  depend  on  the  hardware  implementation. 


4.5.1  Background  and  Quantum  Noise  Limited  Operation 
The  performance  is  approximated  in  Fig.  3.3  as 


PE(i)  <  NT(i)  (K-l)  exp  {-  [Ab(i)  +  Xg(i)  -  4^0)  ]2T(i)}  (18) 

where 
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2-WAY  ZOOM  ACQUISITION  (A).  HETERODYNE  DETECTION.  FM  SIGNAL.  K  =  4.  PE  =  10"4.  NF  =  1 

TOTAL  AND  PARTIAL  SEARCH  TIMES 


STAGE  0 


STAGE  1 


STAGE  4 


INITIAL  BEAMWIDTH 


Fig.  4.2.  Acquisition  with  zooming  receiver,  heterodyne  detection,  no 
frequency  processing:  (c)  optimal  partition  of  T,  2-way  zooming  strategy  (A) 


SEARCH  TIME  (») 


2-WAY  ZOOM  ACQUISITION  (B).  HETERODYNE  DETECTION,  FM  SIGNAL.  K  =  4.  PE  =  ICH,  NF  =  1 

TOTAL  AND  PARTIAL  SEARCH  TIMES 


Fig.  4.2.  Acquisition  with  zooming  receiver,  heterodyne  detection,  no 
frequency  processing:  (e)  optimal  partition  of  T,  2-way  zooming  strategy  (B) 
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1-WAY  ZOOM  ACQUISITION.  HETERODYNE  DETECTION.  FM  SIGNAL.  K  =  4,  PE  =  lO"4.  NF  =  100 

TOTAL  AND  PARTIAL  SEARCH  TIMES 


Fig.  4.3.  Acquisition  with  zooming  receiver,  heterodyne  detection, 
simultaneous  frequency  acquisition:  (a)  optimal  partition  of  T,  1-way  zooming 
strategy. 
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1-WAY  ZOOM  ACQUISITION,  HETERODYNE  DETECTION.  FM  SIGNAL.  K  =  4,  PE  =  10“*,  NF  =  100 

OPTIMUM  PE  FOR  EACH  STAGE 


INITIAL  BEAMWIDTH 


Fig.  4.3.  Acquisition  with  zooming  receiver,  heterodyne  detection, 
simultaneous  frequency  acquisition:  (b)  optimal  parition  of  PE,  1-way  zooming 
strategy. 
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2-WAY  ZOOM  ACQUISITION  (A).  HETERODYNE  DETECTION.  EM  SIGNAL.  K  =  4.  PE  =  10-*.  NF  =  100 
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Fig.  4.3.  Acquisition  with  zooming  receiver,  heterodyne  detection, 
simultaneous  frequency  acquisition:  (c)  optimal  partition  of  T,  2-way  zooming 
strategy  (A). 
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2 -WAY  ZOOM  ACQUISITION  (A).  HETERODYNE  DETECTION.  EM  SIGNAL.  K  =  4,  PE  =  1(H.  NF-100 

OPTIMUM  PE  FOR  EACH  STAGE 


Fig.  4.3.  Acquisition  with  zooming  receiver,  heterodyne  detection, 
simultaneous  frequency  acquisition:  (d)  optimal  partition  of  PE,  2-way  zooming 
strategy  (A). 
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2-WAY  ZOOM  ACQUISITION  (B).  HETERODYNE  DETECTION.  FM  SIGNAL.  K  =  4.  PE  =  10“*.  NF  =  100 
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INITIAL  BEAMWIDTH 


Fig.  4.3.  Acquisition  with  zooming  receiver,  heterodyne  detection, 
simultaneous  frequency  acquisition:  (e)  optimal  partition  of  T,  2-way  zooming 
strategy  (B). 


ERROR  PROBABILITY 


2-WAY  ZOOM  ACQUISITION  (B),  HETERODYNE  DETECTION.  FM  SIGNAL.  K  =  4,  PE  =  10-*.  NF  =  100 

OPTIMUM  PE  FOR  EACH  STAGE 


Fig.  4.3.  Acquisition  with  zooming  receiver,  heterodyne  detection, 
simultaneous  frequency  acquisition:  (f)  optimal  partition  of  PE,  2-way  zooming 
strategy  (B). 
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The  expression  of  f ^( • )  in  (18)  is  in  the  form  of  (IS)  so  that  the  optimum  is 
given  by  (16).  The  optimal  partition  of  PE  and  T,  as  well  as  the  resulting 
total  T  are  shown  for  the  case  K  -  4,  Aj,o  “  103  *n  *.5(a-f)  and 

Abo  "  106  in  Fig.  4.6(a-f),  for  the  three  acquisition  strategies.  For  large 
noise  count  rates.  (Ab0/NT  >  As),  the  first  stage  of  the  search  is 
dominant  so  that  the  three  strategies  are  equivalent.  For  small  noise  counts, 
the  first  stage  is  also  dominant  for  a  two-way  acquisition.  The  optimal 
strategy  for  a  one-way  acquisition  in  low  background  noise  however  is  to  spend 
approximately  the  same  time  at  each  stage  of  the  zooming.  The  two-way 
acquisition  is  therefore  faster  in  this  case  by  a  factor  equal  to  the  number  of 
required  zooming  stages;  see  Fig.  4.7. 

The  background  and  shot  noise  limited  operation  is  difficult  to  reach  in 
practice  without  photomultiplicatlon.  The  photomultiplier  devices  Introduce 
excess  noise  so  that  the  acquisition  times  need  approximately  to  be  multiplied 
by  the  noise  factor  of  the  detector. 

4.5.2  Thermal  Noise  Limited  Operation 

The  expression  for  fi(«)  is  derived  from  Fig.  3.9  as 

Ag(i)  /T(i) 

/  2kt 
/  e2NT(i) 

considering  that  the  thermal  noise  is  dominant.  This  expression  is  in  the  form 
of  (13)  so  that  the  optimality  relationship  is  (14).  Optimal  partition  of  PE 


PE(i)  <  NT(i)(K-l)  erfc 
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Fig.  4.5.  Acquisition  with  zooming  receiver,  background  and  quantum 
limited  direct  detection.  Ago  *  10  •  («>  optimal  partition  of  T,  1-' 
zooming  strategy. 
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Fig.  4.5.  Acquisition  with  zooming  receiver,  background  and  quantum  noise 
limited  direct  detection.  Ago  -  10*:  (c)  optimal  partition  of  T,  2-Way 

zooming  strategy  (A). 
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Fig.  4.5.  Acquisition  with  zooming  receiver,  background  and  quantum  noise 
limited  direct  detection,  Xgo  “  10’:  (e)  optimal  partition  of  T,  2-way 
zooming  strategy  (B). 
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Fig.  4.5.  Acquisition  with  zooming  receiver,  background  and  quantum  noise 
limited  direct  detection,  Ago  53  10  •  (f)  optimal  partition  of  PE,  2-way 

zooming  strategy  (B). 
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Fig.  4.6.  Acquisition  with  zooming  receiver,  background  and  quantum 
limited  direct  detection,  Ago  “  10*:  (•>  optimal  partition  of  T,  1-' 
zooming  strategy. 
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Fig.  4.6.  Acquisition  with  zooming  receiver,  background  and  quantum  noise 
limited  direct  detection,  Ago  -  10®:  (c)  optimal  partition  of  T,  2-vay 

zooming  strategy  (A). 
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2-WAY  ZOOM  ACQUISITION  (A).  DIRECT  DETECTION.  K  =  4,  PE  =  10"4,  A.^  =  10® 
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Fig.  4.6.  Acquisition  with  zooming  receiver,  background  and  quantum  noise 
limited  direct  detection,  XBo  -  10* :  (d)  optimal  partition  of  PE,  2-way 

zooming  strategy  (A). 
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2 -WAY  ZOOM  ACQUISITION  |B),  DIRECT  DETECTION,  K  =  4,  PE  =  10“*  =  108 

OPTIMUM  PE  FOR.  EACH  STAGE 


and  T  are  given  in  Fig.  4.8(a-f )  for  the  case  K  ■  4,  Z  -  10®  for  the  three 
acquisition  strategies.  The  first  stage  of  the  zooming  is  dominant  for  a 
two-way  acquisition  whereas  all  the  stages  have  the  same  importance  for  a 
one-way  acquisition.  The  one-way  acquisition  times  are  therefore  longer  by  a 
factor  equal  to  the  number  of  zooming  stages  as  is  shown  on  Fig.  4.9. 
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Fig.  A. 8.  Acquisition  with  zooming  receiver,  thermal  noise  limited  direct 
detection,  input  impedance  -  106  ft:  (b)  optimal  partition  of  PE,  1-way 
zooming  strategy. 
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Fig.  4.8.  Acquisition  with  zooming  receiver,  thermal  noise  limited  direct 
detection,  input  impedance  =■  106  ft:  (e)  optimal  partition  of  T,  2-way 
zooming  strategy  (B). 
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Fig.  4.9.  Acquisition  times  for  zooming  receiver,  thermal  noise  limited  direct 
detection,  comparison  between  three  zooming  strategies.  Input  impedance  - 
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SUMMARY  AND  CONCLUSIONS 


5.1  Introduction 

Ue  review  In  this  chapter  the  results  we  have  obtained  for  the  spatial 
acquisition  performance  of  a  space  optical  communication  link.  Approximate 
expressions  of  the  acquisition  time  for  different  transmitter  and  receiver 
strategies  are  given  in  Sec.  5.2.  We  suggest  what  the  best  receiver  structure 
and  acquisition  strategy  may  be  in  Sec.  5.3. 
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5.2.1  Parallel  Receiver  Operation 

For  direct  detection,  Eq.  (3.6)  gives  the  perforaance  for  background  noise 
and  quantum  noise  limited  operation  and  leads  to  the  approximation 


T 


(2) 


Photomultiplication  detectors  usually  need  to  be  used  to  reach  the  above  noise 
figures  so  that  excess  noise  has  to  be  taken  into  account. 


T  a 


(3) 


Thermal  noise  limited  operation  results  in  an  acquisition  time  given  by 
inverting  Eq.  (3.9)  and  approximated  by 
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where  leqU  ”  kT/Ze  *  10  / Z  is  a  parameter  which  may  be  compared  to  A^  in 

Eq.  (2).  The  quantum  noise  limit  in  Eq.  (3.7)  gives 
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for  an  error  probability  of  10' 

For  heterodyne  detection,  Eq.  (2.17)  corresponds  to  the  performance  of  the 


acquisition  when  the  signal  has  no  amplitude  fluctuations  and  is  inverted  as 


tcmsaaaan 
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The  receiver  filter  baadwldth  WF  has  the  same  role  as  X^  In  Eq.  (2). 

Equation  (2.21)  corresponds  to  the  unrealizable  perfect  frequency  acquisition 
and  gives 


T  -  t— 
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for  a  random  frequency  f0  and 


(8) 


for  an  a-priorl  known  frequency  fQ. 


5.2.2  Serial  Receiver  Operation 

In  the  case  of  a  sequential  search  of  the  M  spatial  slots  by  K  detectors 
at  the  receiver,  the  figures  given  above  for  parallel  receiver  operation  are 
just  Increased  by  a  factor  of  M/K. 


5.2.3  Zooming  Receiver  Operation 

We  have  seen  in  Chapter  4  that  in  most  of  the  cases,  the  acquisition  time 
is  similar  for  one-way  and  two-way  zooming  and  corresponds  to  the  first  stage 
of  the  search.  The  one-way  zooming  is  much  simpler  to  implement  and  should 
therefore  be  preferred.  When  we  take  into  account  that  each  detector  covers 
M/K  spatial  modes  at  the  first  stage,  we  obtain  the  following  approximations 
for  the  acquisition  time  T. 


Direct  detection,  background  noise  limited  receiver 


T  *  P 


[100  M  ^bo  50 
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-  Direct  detection,  thermal  noise  limited/receiver 


100  equ 


(10) 


-  Heterodyne  detection,  FM  signal 
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For  background  noise  limited  direct  detection,  i.e.,  when  the  first  term 
is  dominant  in  (3)  and  (9),  the  zooming  performance  In  (9)  is  Identical  to  the 
serial  performance  obtained  by  multiplying  (3)  by  M/K.  This  may  also  be  seen 
from  the  comparison  of  Figs.  3.4  and  4.7.  For  heterodyne  detection,  the 
zooming  acquisition  times  given  in  Eq.  (10)  are  worse  than  for  the  serial 
acquisition,  due  to  the  low  efficiency  of  multimode  mixing.  This  Is  obvious 
from  the  comparison  of  Figs.  2.7  and  4.4. 


5.3  Discussion 

The  Issue  of  choosing  the  right  acquisition  strategy  for  a  specific 
situation  depends  on  various  factors  such  as  performance,  hardware  complexity, 
weight,  cost  and  device  availability.  The  choice  may  therefore  be  dependent  on 
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evolving  factors  such  as  technology  state-of-the-art ,  so  that  no  permanent 
answer  exists.  We  will  compare  different  strategies  and  point  out  the  reasons 
why  one  should  be  preferred  to  the  other.  Among  the  various  choices  that  have 
to  be  made,  we  will  address  the  following.  The  illumination  of  the  uncertainty 
zone  may  be  parallel  or  serial,  with  the  illumination  beamwidth  left  as  a 
parameter  for  serial  illumination.  The  receiver  is  usually  difficult  to 
Implement  with  a  parallel  structure,  due  to  hardware  complexity;  the  choice  has 
therefore  to  be  made  between  serial  and  zooming  receiver  structures.  The 
detector  technology  may  be  heterodyne  or  direct  detection,  with  or  without 
photomultiplication. 


5.3.1  Receiver  Technology 

The  comparison  between  the  different  technologies  may  be  simplified  to  the 
following  comparison  for  low  incident  power.  We  approximate  (3),  (4),  (6)  by 
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where  A 
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for  a  PIN  diode  receiver  and  an  amplifier  of  input 
impedance  Z 


A^F  for  an  APD  or  PMT  receiver  with  noise  factor  F,  in  a 

background  noise  limited  operation 


for  heterodyne  detection  with  a  filter  bandwidth  of 
WF 


It  is  clear  that  direct  detection  without  photomultiplication  has  always  the 
worst  performance.  If  Aj,0  ■  106  counts  per  second  and  per  spatial  mode  and 
F  ■  10,  then  with  heterodyne  detection,  the  receiver  needs  to  filter  in 
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bandwidths  of  10  MHz  to  achieve  a  performance  comparable  to  the  performance  for 
photomultiplled  direct  detection.  This  last  technique  has  therefore  the 
advantage  for  achieving  the  same  performance  with  a  much  simpler 
Implementation.  If  progress  Is  done  In  laser  stabilization  and  tuning  however, 
heterodyne  detection  may  gain  the  advantage  over  direct  detection  because  It  Is 
almost  Insensitive  to  background  noise. 

For  larger  Incident  powers,  l.e.,  when  the  second  term  is  dominant  In  (3) 
and  (6),  the  two  technologies  are  comparable  for  F  •  2.  The  optimal  quantum 
limited  performance  given  In  (5)  and  (7)  or  (8)  give  the  advantage  to  direct 
detection. 

Although  heterodyne  detection  can  be  superior  to  direct  detection  for 
signal  demodulation  with  very  large  data  rates,  this  advantage  does  not  apply 
to  spatial  acquisition  when  there  Is  frequency  and  phase  randomness.  Moreover, 
heterodyne  receivers  have  twice  the  noise  bandwidth  of  homodyne  receivers  but 
homodynlng  Is  not  practical  for  spatial  acquisition. 

5.3.2  Receiver  Strategy 

The  optimum  receiver  strategy  is  the  parallel  receiver  in  all  cases,  but 
receiver  arrays  of  104  elements  and  the  accompanying  electronics  are 
unrealistic.  For  heterodyne  detection,  the  sequential  receiver  operation 
outperforms  the  zooming  operation.  For  direct  detection,  the  zooming  and 
serial  operations  give  comparable  performance  for  low  signal  to  noise  ratios, 
when  the  Illumination  beams  are  wide.  Both  receivers  are  then  background  noise 
limited.  When  the  Illumination  beams  are  narrowed  down,  the  serial  receivers 
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reach  the  quantum  limited  performance  earlier  than  the  zooming  receivers;  see 
Fig.  5.1. 


From  expression  (3),  the  serial  receiver  is  background  limited  when 
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From  expression  (9),  the  zooming  receiver  is  background  limited  when 
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Background  noise  limited  operation  is  therefore  obtained  for  larger  values  of 
for  the  zooming  receiver. 

The  expressions  derived  here  therefore  give  the  advantage  to  the  zooming 
receiver  operation.  Caution  should  be  taken  however  because  the  background 
noise  model  Includes  no  spatial  variations.  This  assumption  is  not  verified  in 
practice  and  a  quantum  noise  limited  operation  may  prove  much  more  robust  than 
a  background  noise  limited  operation. 

In  all  cases  except  for  thermal  noise  limited  direct  detection,  the  first 
stage  of  a  zooming  strategy  is  predominant.  Since  the  first  stage  is 
Independent  on  the  zooming  strategy,  the  most  simple  strategy  may  be  chosen. 
One-way  zooming,  which  consists  of  illuminating  the  complete  initial 
uncertainty  zone  during  the  total  search  time,  is  therefore  the  advisable 
zooming  strategy.  Moreover,  the  travel  time  of  light  from  A  to  B  needs  to  be 
added  between  the  stages  in  a  cooperative  two-way  zooming  procedure.  The 
one-way  zooming  strategy  is  therefore  optimal  for  small  T. 
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5.3.3  Illumination  Strategy 

The  performance  is  always  improved  by  illuminating  the  uncertainty  zone 
sequentially  with  a  diffraction  limited  beam.  The  robustness  of  this  strategy 
with  respect  to  relative  motion  may  not  be  adequate  however.  A  performance 
close  to  optimal  may  be  reached  with  a  wider  illumination  beam  for  a  sequential 
receiver  if  the  two  terras  in  expressions  (3)  or  (6)  have  the  same  order  of 
magnitude.  This  implies  that  increasing  above  5X{j0/^3  ^or  direct 
detection  and  2WF/XS  for  heterodyne  detection  does  not  reduce  the 
acquisition  time  significantly  for  a  sequential  receiver  operation. 

5.3.4  Conclusion 

From  the  present  analysis,  it  turns  out  that  two  acquisition  strategies 
are  advisable.  They  both  use  a  receiver  of  photomultiplied  direct  detectors 
and  a  serial  illumination  scheme.  The  receiver  operation  may  be  a  serial  or  a 
one-way  zooming  operation  depending  on  robustness  considerations,  with  the 
advantage  to  zooming  over  serial  for  the  performance. 

Arrays  of  APD  are  not  yet  realizable  because  of  the  large  bias  voltages 
involved.  This  issue  may  be  circumvented  by  linking  the  optical  system  focal 
plane  to  a  number  of  separated  APD  by  an  optical  fiber  bundle.  A  different 
photomultiplied  receiver  array  may  be  obtained  by  combining  a  microchannel 
plate  and  a  photodiode  array. 

Figure  5.1  shows  curves  of  acquisition  times  for  the  two  above  mentioned 
strategies,  for  K  =  10  and  K  =  100.  The  noise  factor  F  is  taken  as  10  and  the 
background  noise  count  rate  per  spatial  mode  Xjj0  as  106  s"1. 
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The  assumption  that  the  spatial  parameter  is  discrete  valued  is  not 

verified  in  practice.  If  the  incoming  signal  is  centered  on  a  corner  of  the 

detector  array,  the  signal  power  is  split  to  four  detectors.  The  signal  power 

has  therefore  to  be  increased  by  a  factor  of  4. 

If  shorter  acquisition  times  are  required,  it  is  possible  to  improve  the 

search  by  providing  a  beacon  laser  of  power  P1^  »  on  one  of  the 

satellites,  say  A.  Satellite  B  will  therefore  acquire  the  position  of  A  with  a 

L  LA 

search  time  given  by  the  above  expressions,  with  P^  replaced  by  PA  .  The 
acquisition  time  decreases  at  least  as  fast  as  the  inverse  of  the  laser  power 
P1^.  The  acquisition  of  the  position  of  B  from  A  is  afterwards  improved 

because  £lP  is  reduced  to  the  diffraction  limit,  so  that  B  may  concentrate  the 
A 

smaller  power  of  its  laser  into  a  smaller  zone.  An  asymmetric  spatial 
acquisition  may  therefore  be  dominated  by  the  fastest  search,  if  the  two 
searches  are  done  successively. 
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APPENDIX  1 


STATISTICS  AFTER  ENERGY  DETECTION  IN  A  BANDWIDTH  WF 


The  Incoming  signal  is  processed  by  the  detector  shown  in  Fig.  Al. 

Signal  Model  (a):  GRP  with  rectangular  PSD,  power  P8 

Bandwidth  W®  one-sided  included  in  WF 
AWGN  of  PSD  N0/2 

Signal  and  noise  statistically  independent. 

The  spectral  density  of  the  signal  x  is  shown  on  Fig.  A2.  The  output  y  of 
the  bandpass  filter  has  a  spectrum  Syy(f)  which  is  also  shown  on  Fig.  A2. 
Furthermore, 


Ryy(0)  ’  L  syy(f)  df  “  pB  +  nomF  “  s V  +  NowF  * 

where  Ss/2  -  PS/2WS  is  the  signal  double  sided  PSD  height. 

The  signal  y  is  squared  and  the  resulting  signal  z  Is  characterized  by 


m 

z 


R 

yy 


(0) 


Rzz(t) 


S 

zz 


(f) 


R2  (0)  +  2  R2  (t) 

yy  yy 

R2  (0)  6(f)  +  2  S  (f)  *  S  (f) 

yy  yy  yy 


The  squared  signal  z  is  Integrated  during  T  seconds  and  the  result  R  is 
characterized  by 


E(R)  -  T  b  -  T  R  (0)  -  P  T  +  N  W  T 

'  z  yy  o 


MR2)  -  7  S„<«  T2  2  Of 


RyytOjT2  +  2  /  [Syy(f)  '  Syy ( t ) )  ^  ^ 


sin  irfT 
irfT  > 


If  Syy(0  is  smooth  over  frequency  intervals  of  the  order  of  magnitude  of 
1/T,  the  convolution  product  may  be  taken  out  of  the  integral 


E(R2)  -  Ryy(0)  T2  +  2  Syy(f)  *  syy(f) 


+,“  J2  /sin  nft\  2 


f*0 


/  T2  (i 

“09  ' 


+0O  2 

Var  (R)  -  E(R2)  -  E2(R)  -  2T  /  [Syy(f)] 


2  /N  s\2 

-  2T  (2/  -  2W8)  +  J  (2W8) 

g2 

-  N2  WFT  +  2N  P®T  +  ■—  T  . 

0  0  w8 

The  random  variable  Ris  a  /  variable  of  order  approximately  2WFT.  If  the 
order  is  larger  than  100,  the  statistics  are  accurately  modeled  by  Gaussian 
statistics. 

In  the  presence  of  the  signal,  we  have  therefore 
R  •  N(ma,  o2) 


A- 4 


a  -  K  W  T  +  P  T 
s  o 


o2  -  N  2WFT  +  2  N  P*T  +  - —  T 
s  o  o  ^ 


When  only  noise  la  present,  the  statistics  are  obtained  from  the  above  by 
inserting  P8  -  0: 


R  -  N  (a  ,  o  ) 
n  n 


d  -  N  W  T 

n  o 


J2  -  N  2WFT 


Signal  Model  (b): 


Signal  with  frequency  noise  but  no  amplitude  noise 
Bandwidth  W®  one-sided  included  in 
AWGN  of  PSD  M0/2 

Signal  and  noise  statistically  Independent. 


Ue  conserve  the  notation  used  for  signal  model  (a). 


x  -  characterization 


x(t)  -  8 ( t )  +  n(t) 


Elsz(t)J  -  P8 


y  -  characterization 


y(t)  *  s(t)  +  n^(t) 


where  nf(t)  is  the  filtered  version  of  n(t),  and  assuming  the  signal  is 
completely  in  the  filter  bandwidth. 


z  -  characterization 


z(t)  -  s2(t)  +  a(t)  nf(t)  +  n2(t) 

»z  -  E  [ a2 ( t ) ]  +  2  E[s(t) ]  E[nf (t)]  +  E[n2(t)] 

-  P8  +  N  WF 
o 


R  -  characterization 


-  /  mz(t)  dt  -  P8T  +  N  W^T 

o 

T  T 

Var  (R)  ■  Var  /  2s(t)  n^(t)  dt  +  Var  /  n^(t)  dt 

o  o 

The  first  tern  is  the  variance  of  a  zero-mean  variable  and  therefore  equal  to 

2 

its  mean  square.  An  expression  for  the  second  tern  is  given  by  oq  in  model  (a) 
above. 


Var  (R) 


N^T  +  4  /  dt 

o  ' 
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T 

/ 


o 


du  E[s(t)s(u) 1  E[nf (t)nf (u)J 


2  F 

■  irw  t  + 
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T 

I 

o 


dt 


T  N 

/  du  Rgs(t  -  u)  6(t  -  u) 


N2WFT  +  2N  P*T 


In  this  case  also,  R  is  the  sum  of  approximately  2Wr T  independent  random 


variables.  If  Wft  >  100,  we  have: 

R  -  N(m  ,  of) 
n  n 

m  -  PST  +  N  WFT 
n  o 

o2  -  N^T  +  2N  PST 
no  o 

Signal  Model  (c):  Random  phase  sinusoid 

AWGN  of  spectral  height  NQ/2 
Signal  and  noise  statistically  Independent. 
We  conserve  the  notation  introduced  for  signal  (a) 


*  (t)  -  cos  2  irf  T  +  n(t) 

o 

y  (t)  -  SIP  cos  2nf  t  +  n.(t) 

or 

z  (t)  ■  y2(t)  -  2P8  cos2(2irf  t)  +  2  /  2P8  cos(2irf  t)  n  (t)  +  n?(t) 

o  o  t  t 

T 

^  -  /  m^Ct)  dt  -  P®T  +  NqWFT 

o 


Var  [ R]  ■  Var  /  2  /  2P8  cos(2nf Qt)nf(t)dt  +  Var  /  n2(t)dt 


From  the  projection  of  the  flat  spectral  density  noise  on  one  coordinate  of  an 
orthogonal  expansion,  we  know  that 


Var 


/  ^  cos(2*fot)  n(t)dt 


The  first  term  in  the  expression  of  Var  [R]  is  therefore  4PST  NQ/2.  The 
2 

second  term  is  from  signal  model  (a)  and  therefore. 


c £  -  N^W^T  +  2N  PST 


APPENDIX  2 


The  LRT  is  obtained  as: 


(R,  -  R0)<R,  +  Ro  -  2RT)  ^  0 


or 


H. 


Is  -  rt|  i  K  -  si 


where 


R_  - 


2  2 
m  a  -mo 
ns  s  n 

2  2 

a  -  a 
s  n 


(A2-3) 


(A2-4) 


(A2-5) 


The  resultant  partition  of  the  decision  plane  (Rg,  Rj)  into  (ZQ,  Zj)  is 
shown  in  Fig.  A3. 

The  implementation  of  the  optimal  rule  (A2-4)  implies  the  knowledge  of 
R.p  related  to  the  signal  parameters. 

The  intuitive  decision  rule 


H 

R1  * 
H 


1 


(A2-6) 


is  only  suboptlmal,  but  is  close  to  (A2-4)  when  R^  1#  small.  The  range  of 

Rt  is  from  -<*>  when  og  ■  on  to  when  og  >  Ojj.  The  signals 

considered  In  Chapter  3  result  in  the  following  values.  For  Signal  Model  b) , 

rt  *  ojjj/2.  For  Model  a),  mn/2  <  RT  <  n^.  We  see  therefore  that 

the  optimal  decision  rule  is  also  dependent  on  the  signal  model.  We  therefore 

consider  only  the  suboptlmal  rule  given  in  (A2-6). 


The  sufficient  statistic  for  the  suboptiaal  rule  in  (A2-6)  is 


H. 
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i  *  N(m  -  m  ,  a2  +  o2)  Ih 
s  n*  s  n  I  1 

l  *  N(-m  +  m  ,  a2  +  a2)  Ih 
s  ns  n  I  c 
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M-ary  Signaling  The  optimum  decision  rule  is  to  find  the  largest 
|Ri  -  Rqj  where  RT  is  given  in  (A2-5).  The  suboptiaal  rule  of  choosing 
the  largest  %  has  an  error  probability  which  is  upper  bounded  by 


PE  <  (M  -  1)  erfc 


m  - 
8 


n 


o2  + 
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APPENDIX  3 

COMPOSITE  M-ARY  DETECTION  WITH  AN  N-ARY  UNWANTED  PARAMETER 
APPLIED  TO  ORTHOGONAL  RANDOM  PHASE  SINUSOIDS 


One  out  of  N  x  M  orthogonal  signals  may  be  received  with  equal  a  priori 
probabilities.  This  set  of  possible  signals  is  divided  into  M  classes  of  N 
signals.  The  relevant  problem  is  to  detect  the  class  of  the  received  signal. 
The  unwanted  parameter  specifies  the  particular  signal  in  the  class. 
Mathematical  formulation: 

Under  hypothesis  H^ j ,  £( t)  ■  S^j(t)  +  n(t)  (A3-1) 

All  SjjCt)  are  orthogonal,  ii(t)  is  AWGN. 

The  index  i  specifies  the  relevant  parameter  and  j  the  unwanted  parameter.  For 

our  particular  problem,  i  specifies  the  spatial  position  and  j  specifies  the 

frequency  and  time  of  arrival  of  the  signal.  The  analysis  presented  below  is 

kl 

more  general  however.  If  C^  ■  cost  of  declaring  hypothesis  kl  when  signal  i^ 
is  present,  then  the  cost  assignment 


(A3-2) 


emphasizes  no  cost  is  incurred  for  j  *  1  as  long  as  i  ■  k.  The  optimal 
decision  rule  is  to  find 


(A3-3) 


A- 13 


and  declare  hypothesis  Hp  If  the  maximum  is  attained  for  1  ■  p.  In  the 
particular  case  where  Sjj(t)  are  orthogonal  random  phase  sinusoids,  the 


observation  Is  reduced  to 


fijc  + 


(A3-4) 


where  Rj_js  and  Rijc  represent  the  sine  and  cosine  Integrated  signal  St j . 


The  conditional  statistics  are 


-xIh^IV  • 


4  "".'".(il  a-l  I^"  +  E)1  " 

■5“  e  I  I « —  n  n 


(v  »i  X|ai  (u-5) 


The  optimal  decision  rule  is  to  find: 


r  T  /2SE  „  \ 
nax  >  I  I  s —  X,  4 1 
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(A3-6) 


And  declare  hypothesis  Hp  if  the  maximum  is  attained  for  i  «  p.  For  large 
SNR,  IQ(x)  has  an  exponential  dependence  on  x.  The  suboptimal  rule 


max  max  X 

i  j  5 


(A3- 7 ) 


is  therefore  close  to  optimal  for  large  SNR,  Furthermore,  (A3-7)  is 
independent  of  the  signal  parameters. 


Exact  error  probability  for  the  suboptimal  receiver  (A3 -7 ) 


PE  -  1  -  PAC  -  1  -  PAC 


(symmetric  problem) 
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P  +  (N  -  l)  P. 
a  d 


We  now  evaluate  Pa  and  P^. 
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(A3-8) 


This  is  the  exact  error  fornulation  for  the  suboptiaal  receiver  in  (A3-7).  The 
sum  in  (A3-8)  is  upper  bounded  by  its  first  term. 


PE 


(A3- 9) 


The  formula  (A3-9)  is  a  union-type  bound;  an  error  may  occur  only  when  the 
received  signal  projected  on  one  of  the  (N  -  1)  non-signal  classes  exceeds  the 
projection  on  the  transmitted  signal. 


APPENDIX  4 

BOUNDS  TO  POISSON  DETECTION  ERROR  PROBABILITIES 


a)  Binary  Detection 

The  two  discrete  valued  independent  random  variables  r^  and  rQ  are 
conditionally  Poisson  distributed  with  means  m ^  and  m0  on  Hlt  or  m0  and 
on  Hq(  with  >  mQ.  Generally  mQ  -  nB  is  the  expected  number  of 
background  counts  and  idj  *  nB  +  ng  is  the  expected  number  of  background  + 
signal  counts.  The  optimal  decision  rule  for  MEP  and  equal  a  priori 
probabilities  is 
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The  error  probability  is  given  by 
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(A4 


assuming  that  equal  counts  results  in  an  error.  Equation  (A4-1)  may  also  be 
expressed  as  [7]. 


PE 


I  (2/m_m1 ) 

p  0  1 
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b)  M-ary  Detection 

The  optimal  decision  rule  is  again  the  largest  count  selection.  The  error 
probability  is  given  by 


PE 


1/  \  n 
e  (mp 

n! 


n-1 

M-0 
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(A4-8) 


The  M-ary  error  probability  is  always  upper  limited  by  the  union  bound,  using 
any  binary  error  bound.  Considering  (A4-6)  for  the  binary  error,  we  have: 


PE  <  (M  -  1)  e 


(A4-9) 


This  bound  is  tight  when  the  (M  -  1)  possible  binary  errors  are  lowly 
correlated.  In  the  case  of  low  however,  the  binary  errors  are  highly 
correlated  and  (A9-9)  may  over-estimate  PE  by  a  factor  of  up  to  (M  -  1). 

Another  M-ary  error  bound  is  found  by  noting  that  no  error  occurs  if  the  signal 
count  is  higher  than  the  sum  of  all  noise  counts. 
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This  bound  gives  the  exact  PE  for  *  0 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  FADE  (Timm  Data  Flir.e 


REPORT  DOCUMENTATIOM  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

f 

1 

Z.MVTACCEMMRR0. 

3.  RECMERTS  CATAUW  RUM0EA 

ESD-TR-84-016 

4.  TITU  (end  Subtitle) 

1.  TYPE  OF  REPORT  SPEIR00  COVERED 

Spatial  Acquisition  Algorithms  and  Systems  for  Intersatellite 

Technical  Report 

Optical  Communication  Links 

A  PERFORMED  PRO.  REPORT  RDHOER 

Technical  Report  667 

T.AIfTNOAM 

1.  CONTRACT  OR  RRART  RUMRERfV 

Patrick  Van  Hove  and  Vincent  W.S.  Chan 

F19628-85-C-0002 

t.  FEAF0AMM6  OASARUATKM  RAME  ARO  ADDRESS 

11.  PR00RAHI  EtERKRT.  PROJECT.  TAM 

Lincoln  Laboratory,  M.I.T. 

AREA  R  WORE  INHT  HUMKRt 

P.O.  Box  73 

Program  Element  No.  6343 IF 

Lexington,  MA  02173-0073 

11.  C0RTR0UIR0  OFFICE  RAISE  ARO  ADORE  IS 

11.  REPORT  DATE 

Air  Force  Systems  Command,  USAF 

Andrews  AFB 

27  November  1984 

IS RUMRER  OF  PARES 

|  Washington,  DC  20331 

132 

1  14.  MORITOAIHG  AOEHCY  RAME  A  ADOAESS  (if  different  from  Controlling  Office) 

It.  SECUIRTY  CLASS,  (ofthu  report) 

Electronic  Systems  Division 

Unclassified 

Hanscom  AFB,  MA  01731 

Ik.  DECUSMFICATWR  00WR6RADNM  SCHEDULE 

II.  OISTMBUTHM  STATEMEHT  (of  tbit  Report) 


Approved  for  public  release;  distribution  unlimited. 


IT.  DISTMSimOD  STATEMERT  (of  the  abtrart  entered  in  Black  SO,  if  different  from  Report) 


IS  SUFFUNESTMT  ROTES 

None 


IS  KEY  WORDS  (Continue  on  reveree  *1 do  If  nummary  and  identify  Ay  Mack  number) 

optical  communication  4  laser  communication  Cm-  „ 

spatial  acquisition  '  >  satellite  cross  links 

_ 

2t.  ASSTRACT  (Continue  on  mono  tide  if  nummary  and  identify  by  block  number) 

\ 

\ 

Spatial  acquisition  strategies  and  technologies  are  analyzed  for  optical  communication 
systems.  Theoretical  analysis  is  followed  by  examples  and  comparisons  to  aid  system  design. 
Examples  are  given  for  state-of-the-art  technology  appropriate  for  space  applications. 

(  1 1  C\  t  i'  O  v  -  .  4>  is- J  Co  - 


DO  fowl  )473  toiTIOi  OF  1  MOV  H  IS  OMOUTE  UNCLASSIFIED 

1  Jsa  73 


MCURITY  ClASSIFICATION  OF  THU  FADE  (When  Dm*  famR 


END 

FILMED 

4-85 

DTIC 


